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Abstract 

In this paper, inspired by the holographic dual of the entanglement entropy, 
we consider the behaviour of extremal, codimension two, spacelike surfaces in 
the background of three and four dimensional charged boson stars in asymp- 
totically anti-de Sitter spacetime. We find conditions for which families of 
minimal area surfaces fail to contain the entire bulk spacetime and construct 
a phase diagram showcasing the transition between regimes. In addition, we 
use the relation between the star's mass and the central density of the scalar 
field to argue for a possible instability of such hollow solutions. Finally, we 
discuss the consequences of our findings for the study of holographic duals of 
reduced density matrices. 
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1 Introduction 



In this paper we construct solutions of asymptotically AdS boson stars coupled to 
a U{1) gauge field in 3 and 4 dimensions, compute the star's charge and mass as 
functions of the scalar field central density and study the behaviour of extremal 
surfaces in these backgrounds. In particular, we determine conditions under which 
minimal area, codimension 2 spacelike surfaces fail to cover the whole spacetime. To 
understand the motivation for doing so, let us digress for a moment and consider 
some intricate open questions regarding the holographic principle in anti de Sitter 
spacetimes. 

The holographic principle is a remarkable idea relating theories with gravity to 
theories without gravity [U [2] . The best understood example of holography is the 
AdS/CFT correspondence which proposes a one to one correspondence between a 
gravitational theory on anti-de Sitter space and a strongly coupled, large N, confor- 
mal field theory living on the conformal boundary of AdS [3ii4j. This correspondence 
has been extensively used to study a variety of strongly coupled field theories and has 
led to the formulation of successful formalisms such as AdS /Q CD and AdS/CMT[^, 
however, the potential use of the gauge gravity duality to tackle problems in quantum 
gravity has yet to show comparable progress. 

The intrinsic non local aspect of the bulk theory is particularly challenging when 
trying to understand the exact connection between bulk and boundary degrees of 
freedom. It stands as one of the major barriers in the way of the gauge / gravity 
duality fulfilling its potential of answering long standing problems in quantum gravity 
by recasting them in field-theoretic language. In order to improve this situation, the 
issue of understanding the map between bulk and boundary degrees of freedom must 
be addressed. Optimistically, we should be able to express any variation of the 
boundary theory state in terms of a well defined perturbation of the bulk geometry; 
conversely, any changes in the bulk should be related to a particular perturbation of 
the boundary state in a one to one fashion. 

An important question related to establishing the dictionary between bulk and 
boundary is to understand what information about the bulk is contained in a certain 
region of the boundary field theory. In other words, let a local field theory defined 
on the conformal boundary B of a spacetime be in a generic state p. Assume this 
theory has a well defined gravity dual and consider, with respect to some arbitrary 
time slicing, a spacelike subregion A of the boundary. Which portion of the bulk 
R{A) is dual to the reduced density matrix pjf], is a question closely related to the 

^For an extensive review, refer to [SI [5] ■ 

^Where pA is the partial trace of p over the complement of A, A. 
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issue of characterizing the bulk-boundary degree of freedom map. 

The Ryu-Takayanagi proposal [7j suggests that the entanglement entropy of a sub- 
region y4 of a field theory with a well behaved, static, gravity dual, is proportional to 
the area of the minimal area surface anchored at the boundary 6 A of A. Using this 
proposal we can construct a family of minimal surfaces by considering the surface 
dual to the entanglement entropy of pA and all the surfaces dual to the partial traces 
of Pa with respect to arbitrary subregions within A. This family of surfaces defines a 
region in the bulk we will call w(A) and is a possible candidate for R{A) as proposed 
in and further explored in [9p. 

Despite meeting several conditions for a suitable candidate for R{A) [8l|9], when 
A = B, the whole boundary slice, w{B) does not cover the entire spacetime in 
general. It was pointed out in [8J that there are spacetimes for which no minimal 
surfaces reach certain areas of the bulk that are, nevertheless, causally connected to 
the boundary, leading us to the conclusion that, at least for such cases, w{A) cannot 
be the elusive region R{A). 

An explicit example of minimal surfaces that do not cover the bulk is found in 
asymptotically AdS boson star backgrounds. Boson stars are solutions to Einstein's 
equation coupled to a complex scalar field that have attracted the interest of physi- 
cists for over 40 years. Beginning with the work of Kaup [10] and others [11], the 
fields of general relativity, cosmology, and even particle physics have shown great 
interest in fully understanding these solitonic objects. Some of the standard reviews 
are [HI US [H HH] . 

In this paper we compute the physical charge and mass of asymptotically AdS, 
charged boson stars as a function of the scalar field central density, investigate nu- 
merically the behaviour of extremal surfaces on these backgrounds, and determine 
the conditions for which w{A) fails to fully cover the bulk, which we will call hollow 
phases. Furthermore we compare the conditions for the existence of hollow phases 
with that of known physical instabilities in four dimensional charged boson star sys- 
tems, and argue that such hollow configurations are likely unstable and therefore 
physically forbidden, thus providing further evidence that w{A) may be the correct 
proposal for R{A). 

This paper is organized as follows: in section 2 we present the action, the metric 
and fields ansatz and equations of motion that follow, in addition we discuss boundary 
conditions and the free parameters we are left with once these are imposed. In section 
3 we explore the relation between the star mass, charge, the scalar field central 
density, and the stability of charged boson stars. We present some numerical results, 

^Formally we should extend the discussion to Da, the causal development of A, however, for 
clarity and objectivity's sake we choose to omit it in the Introduction. 
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compare them to what is known from the hterature, and argue that for a certain 
range of parameters the solutions we find are unstable in D > 4 dimensions. In 
section 4 we review and extend the above discussion as well as outline the numerical 
strategy to extract information about extremal surfaces and present the reader with 
preliminary results. In section 5 we construct phase diagrams displaying the relation 
between the free parameters, the transition between solid and hollow regimes, and 
the transition between stable and unstable configurations. We conclude with a few 
final remarks and future directions in section 6. 



2 Equations of Motion 



We start by considering the Einstein-Maxwell action with a negative cosmological 
constant minimally coupled to a massive complex scalar field in D spacetime dimen- 
sion^ 



(2 

We want to restrict our attention to stationary, spherically symmetric configura- 
tions and allow for electric charges only, therefore we will let the metric be of the 
form 

ds^ = -f{r)de + ^ + r^dnD-2. (2.2) 

and adopt the following ansatz for the scalar and gauge fields 

m = i;{r)e-"'\ (2.3) 
Ao = A{r), (2.4) 
Ai^o = 0. (2.5) 

To include the gravitational contribution from both the scalar and gauge fields 
we write the total energy momentum tensor as 

T,, = Tff + T^f , (2.6) 

the sum of the scalar's energy momentum tensor 

TT^^ = Vi^*Vj^ + Vi^Vj^* - Qij (|V^|^ + , (2.7) 



^ Since the main goal of this paper is to discuss numerical solutions of this action, the Hawking- 
Gibbons surface term can be neglected for it does not alter the equations of motion. 
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where Vj is the covariant derivative (di—iqAi), and the gauge field energy momentum 
tensor 

= F,,F^' - \g,,F\ (2.8) 

Additionally, note that the action [2?T] is invariant under the global U{1) rotation 
\Ef — )■ e*"^', implying the existence of the conserved current 

= tg^\{Vt,mym - (v,.^)^*), (2.9) 

which will act as a source for the gauge field. 

With the above definitions in hand we obtain four linearly independent equations 
of motion from the tt and rr components of Einstein's equation 

^ij ~ 2^ij^ ~ ~ ^ ^~9ij = Tij, (2.10) 

the Klein- Gordon equation 

V^0-m20 = O, (2.11) 

and Maxwell's equations 

V^F'^ = q.P. (2.12) 

In terms of the functions defined in our ansatz, the D dimensional equations of 
motion from the tt and rr components of Einstein's equation are 

- gVA(r)2^/>(r)2 - 2qr^uA{r)^{rY - ^{D - 2)r^-^f{r)g'{r)- 

^{D -3)iD- 2)f{r)g{r) + i(D - 2)(D - l)r^f{r) - r^g{r)A{rf+ 

^{D -3){D -2)f{r) -r^f{r)g{r)'iP'{ry -m\^f{r)'iP{r)^ -r^u^-iPir)^ = 0, 

(2.13) 

and 

gV^A(r)2^/'(r)2 + 2qr^uA{r)ij{rf - ^{D - 2)r^-^f'{r)g{r)- 

^{D-3){D- 2)f{r)g{r) + ]^{D - 2)(D - l)r^f{r) - r^g{r)A{rf+ 

^{D -3){D -2)f{r) +r^f{r)g{r)f{rf -m^r^f{r)^{rf + r^u;^^{ry = 0. 

(2.14) 
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While from the Klein-Gordon equation we find 



2q^rA{rfip{r) + 2{D - 2)f{r)g{r)ip\r) + rg{r)f'{r)ij'{r) + 
rf{r)g'{r)ip'{r) + 2rf{r)g{r)ip"{r) — 2w?rf{r-)ip{r) + 2ruj'^ip{r) = 0, 

(2.15) 

and Maxwell's equations give 

- 2rf{r)g{r)A"{r) - 2{D - 2)f{r)g{r)A'{r) + rg{r)A'{r)f{r)- 

rf{r)A'{r)g'{r)+Aq'^rA{r)f{r)ilj{ry + Aqrujf{r)ilj{ry = 0. 

(2.16) 

These comprise a set of two first order and two second order ordinary differential 
equations and allow for a six parameter family of solutions, one of which is empty 
AdS. However, we are not interested in any generic solution of the above set of 
equations, but only in those that are asymptotically AdS and regular at r = 0. 

By requiring the solution to be regular at the origin, the near r = analysis of 
the equations of motion leads to the following constraints 

g{0) = l, g'iO) = 0, /(0) = /o, f (0) = 0, (2.17) 

A(0) = Ao, A'(0) = 0, V^(0)=V^o, V^'(0) = 0, (2.18) 

leaving three undetermined parameters. We can use the following symmetry of the 
equations of motion 

/^7V, A^^A, co^^co, t^-t, (2.19) 
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to fix A(0) = 1. In addition, the value of /o is chosen to ensure that our solutions 
asymptote to global AdS, that is as r — )■ oo, /(r) — )■ 1 + r^, and similarly for g{r). 
This leaves us with one free parameter, ipo, and the task of studying a one parameter 
family of solutions of equations 12. 13112. 161 

However, before we move forward, we are still left with the issue of guaranteeing 
that the scalar field showcases the correct asymptotic fall off. In general, the large r 
behaviour of ?/'(r) is given by 



with 

\± = \ [{D - 1) ± ^{D-lf + Am^) , (2.21) 

where is the scalar field mass in the lagrangian I2.H and ipi and V'2 are constants. 
For > 0, ilj{r) has a non normalizable term that can be set to zero by picking a 
solution for which tpi = 0, while for (D — 1)^/4 < < both terms are normalizable 
and therefore allowed in principle. 

To pick the desired fall off of V^(r) we will look for the lowest value of u, the phase 
of "^{t), for which we observe ipi = 0. We do so by using a; as a shooting parameter 
and imposing the boundary condition ipi < 10~^°. Once u is fixed we are left with 
only and q (the scalar charge) as free theory parameters. Henceforth in this 
paper we will numerically study the one parameter family of solutions of equations 
l2.13H2.16l their dependence on the solution's parameter tp{0), the central density of 
the scalar field, as well as on the theory's parameters and q. 



3 Mass, Charge and Scalar Central Density 

Before we start the discussion of how the solutions of equations 12. 13112. 161 and the 
families of minimal surfaces depend of the free parameters discussed above, lets take 
some time to look at how the total mass and total charge of the star are calculated 
and how they depend on the central density of the scalar field, ^^(O). 

The metric (12.21) describes a spherically symmetric body whose mass can be 
extracted from the asymptotic behaviour of the metric function g{r). At large r we 
expect the metric to approach that of a charged, massive star in AdS with no scalar 
field, i.e.: 

2M 

g{r) >l + r - + (3.1) 

where M is the mass of the star and Q its charge, which is simply the total number 
of scalar particles times the charge q, i.e.: 

Q = q [ d^-^xJ^^. (3.2) 



However, the existence of a scalar field allows for the possibility of scalar hair, 
therefore we should expect that our metric, while still approaching 13 . 1 1 at large r, will 
receive higher order corrections. Ergo we choose to write the general metric as 



2 2M(r) Q{rf 



9{r) = l + r + (3.3) 
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where the large r behaviour of the function M(r) is given by 

M(r) ~ M + C(l/r"), (3.4) 

with a being is a positive constant, and similarly for Qir). 

The behaviour of the mass and charge of a boson star as a function of the scalar's 
central density has been extensively studied in the literature both for asymptotically 
fiat and asymptotically AdS spacetimes [161 IIZ]- particular, in AdS space, the 
existence of a variety of new solutions and the presence of a zero temperature phase 
transition have recently been shown to exist (TTJ [TS] . Furthermore the stability of 
boson stars has also been the subject of numerous studies and found to correlate with 
certain aspects of the behaviour between mass, charge and central density [TBI 112] 
as we discuss below. 

Since later in this paper the question of whether the solutions we find displaying 
hollow configurations of minimal area surfaces are physically stable will be particu- 
larly important, we should take a closer look at the behaviour of the mass and charge 
of our solutions as a function of the central density '?/'(0) and compare it to what is 
known from the literature. Our analysis here will rely solely on previously known 
results and educated guesses, so we will refrain from a formal study of the stability 
of the solutions as it lies outside the scope of this work. 

For D > 4 dimensions it has been found that for models with no gauge field in 
both zero and negative cosmological constants background, and with a gauge field 
in a zero cosmological constant background, the star mass as a function of the scalar 
field central density M(?/^(0)) reaches a maximum value for a finite '?/'(0) = ^/'c- Even 
though the numerical value for the maximum mass and ^/'c change for each case, they 
correspond to the threshold between dynamically stable ('?/'(0) < t/'c) and unstable 
('?/'(0) > ■j/'c) regime^. 

The D = 3 case has attracted considerably less attention and consequently, to 
the extent of the author's knowledge, no formal result is available. Nevertheless, in 
the context of strongly self interacting boson stars, there is evidence for the existence 
of a maximum mass as well as some discussion regarding the positive binding energy 
of these objects possibly being an indicator of instabilities [201 EI]- Despite the 
existence of such partial results we will refrain from making any statements regarding 
the stability of the three dimensional solutions. 

With the above discussion in mind we numerically solve equations 12. 13112. 16| 
compute the physical mass and charge of the configurations using 13.41 and 13.21 for 

®At maximum mass (^/'(O) = V'c) the pulsation equation arising from the analysis of the time 
evolution of infinitesimal radial perturbations has a zero mode indicating that il)^ is a boundary 
between stable and unstable equilibrium configurations [19) . 
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Figure 1: Plot of the star's mass (orange) and charge (green) versus the central value 
of the scalar field ipQ with = and q = 0.2 in D = A dimensions. 

different values of ip{0) and compare to what is known from the literature. We find 
qualitatively similar results in D = 4 to that of pure boson stars (without gauge 
field), more importantly we observe the existence of a maximum value for M as can 
be seen in figure [1] 

Given the proximity of the models considered in the literature and ours, together 
with the qualitative agreement between results, we are led to believe that our model 
is also unstable past the central density threshold ipc for D > 4 dimensions. We 
will assume this conjecture is indeed true in the reminder of this paper and use it to 
understand better the conditions for the existence of a hollow w{Da) in the upcoming 
sections. 

In D = 3 dimensions we again observe qualitatively similar results to those found 
in [16] (figure [2]). Most notably we would like to higlilight the absence of a maximum 
mass within the central density range investigatecu as it showcases a considerably 
different behaviour than that expected for D > 4 dimensions. 

''Our numerical solution is untrustworthy past the ip{0) = 0.8 mark in 13 = 3 dimensions, 
therefore precluding us from investigating the mass versus central density relation any further. 
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Figure 2: Plot of the star's mass (orange) and charge (green) versus the central value 
of the scalar field ipo with = and q = 0.2 in D = 3 dimensions. 

4 Overlapping Extremal Surfaces 

In the Introduction we raised the question of what portion R{A) of the bulk spacetime 
M. is dual to the reduced density matrix pA for the causal development Da of a given 
region A of the boundary theory. A natural starting point would be to consider the 
causal wedge z{Da) = J^(-D^)n J~(i5^)|^, the region in the bulk a boundary observer 
restricted to Da is causally connected to. Intuitively speaking, we expect z{Da) to 
be at least the minimal portion of the bulk accessible to a boundary observer living 
in Da- The reason being that its complete causal connection to Da allows for the 
detection of any perturbation of the bulk metric inside z{Da)- 

While z{Da) seems to be a natural guess, and it certainly imposes a lower bound 
on the size of R{A), examples for which z{Da) C w{Da) while z{Da) ^ w(Da), are 
easy to finci^ and demonstrate how the region z(Da) alone cannot, in general, hold 
all the information stored in Da- In contrast, for some cases, we can argue that all 
the information in Da must lie within w{Da), the domain of dependence of w{Da), 
for anything outside this bulk region will causally interact with the complement of 

^Where J~^{Da), defined as all points accessible by causal curves arising from Da, is the causal 
future of Da, while J^{Da), defined as the set of points from which Da can be reached following 
a causal curve, is its causal past. 

^See [5] for a thorough discussion. 
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Da, Da- This could lead us to naively expect that w{Da) imposes an upper bound 
on R{A), however, one can construct explicit examples for which this might not be 
true, as we will see below. 

Say that we let the region Da cover the entire boundary, i.e.: Da = B, an 
observer within Da will have access to a full Cauchy surface and, consequently, 
information about the entire past and future of the boundary theory. Since we are 
considering a field theory with a well defined gravity dual, information about the 
past and future of the boundary should extend to information about the past and 
future of the bulk theory as well. In other words, an observer in B that has access 
to all possible boundary physical observables should be able to fully reconstruct the 
dual bulk metric0. 

If we restrict the access of this observer to knowledge of the entanglement entropy 
for any arbitrary region within B, we can ask how much of the bulk he or she can 
probe and, better yet, infer the geometry of [8],[9ll27j. If the bulk theory has a horizon, 
say a spherically symmetric black hole at the origin, it was shown in [28] that, while 
no extremal surface of any co-dimension (or causal lines, for that matter) can probe 
inside the horizon, given a fixed boundary region, co-dimension 2 surfaces probe the 
bulk deeper than higher co-dimension surfaces or causal lines. Now, if instead of a 
black hole at the origin we consider a boson star for example, we can ask the same 
question again: how much of the bulk spacetime can a boundary observer probe with 
extremal surfaces only? Note that now the entire bulk is causally connected to the 
boundary, so we know that z{^) D M., therefore, if w(-Db) fails to fully cover the 
bulk, we will have an explicit example for which z D w while w ^ z. 

Given the above discussion, our goal is to search within the space of solutions of 
charged boson stars in asymptotically AdS spacetime for configurations for which we 
observe overlapping of extremal surfaces leading to a hollow w{Da), therefore ad- 
dressing the question of whether z(Da) 3 w{Da) while w{Da) ^ z(Da) is feasibl^. 

Our setup is both static and asymptotically globally AdS, therefore our boundary 
is a sphere. We will look for surfaces that extend in all polar angles, so, to our 
applications, it suffice to describe them as curves 9{r), with anchor points 6{oo) = 
±6*0 corresponding to the azimuthal boundary coordinates of the start and end points 
(see figures [3l H] and [5]) . 

^"This is demonstrably true for the case of empty AdS [22j [23l [24l [25]. Some interesting re- 
cent discussion highlights the necessity of including non local boundary operator in the boundary 
observer's toolbox [26 

^^Note that failure to find solutions obeying z{Da) D w{Da) while w{Da) ^ ziDA) does not 
indicate this particular phenomena is not possible, in the same way that the mere mathematical 
existence of such configurations is not enough to undermine the candidacy of 'w{Da) to the position 
of R{A) for these could not be physically preferred. 
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Figure 3: A plot of the Penrose diagram of a time slice of multiple extremal sur- 
faces on an asymptotically AdS charged boson boson star background in global co- 
ordinates for = 0, g = 0.1, and ip{0) = 0.2. From top to bottom we have 
6 = ±0.25l7r, ±0.3557r, ±0.4467r, iO.Svr. In this particular case the central density of 
the scalar field is below the threshold iph, therefore there are no degenerate extremal 
surfaces (see figure [6]), in other words, we observe a solid w{Da)- 

To determine whether the entangling surfaces reach the deepest regions in the 
bulk we will analyze the behaviour of r^i,^{6n{r)) for multiple surfaces with distinct 
boundary anchor points 6n{oo) = ±^„, that is, the minimum value of r reached by a 
given entangling surface ^„(r) with fixed boundary anchor points 6n- 

We know that 9{oo)r^^^.^^ (the anchor points of a surface that has rmin as its deepest 
point) covers all values of 6 as we vary rmin from zero to infinity, as a result, for every 
< To < cxD there is at least one extremal surface that obeys r^i^{9{r)) = tq (figure 

ED- 

However, the situation is much different if we focus on surfaces anchored at the 
same boundary points. In this case we expect that only one minimal area surface 
can be anchored at each given 6i (figure l6a|) . While this is true in general (with the 
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Figure 4: Again, a plot of the Penrose diagram of a time slice of multiple extremal 
surfaces on an asymptotically AdS charged boson boson star background in global 
coordinates for = 0, g = 0.1, and ip{0) = 1.2. However, in this example, the 
scalar central density is above the threshold tph and we observe the existence of 
degenerate extremal surfaces for a range of boundary anchor points 6 (see figure [6]). 
In particular, for anchor points 6 = ±7r/2 there are three solutions two of which (blue 
line) lie on top of each other, have minimal area and do not penetrate the dashed 
small circle, while the third (red line) corresponds to a non minimal area extremal 
surface. As discussed in this paper, no minimal area surface penetrates the deepest 
bulk points within the small dashed circle. 

exception of = 7r/2), this need not be the case if we consider extremal rather than 
minimal area surfaces (figure l6b|) . therefore we should not be alarmed if for certain 
boundary anchor points we find that there are multiple distinct extremum surfaces 
anchored to it, in fact, we are most interested in determining when such a phenomena 
happens. To do so we will numerically investigate conditions under which there exist 
multiple extremal surfaces anchored at the same boundary point that do not share 
the same deepest point in the bulk. 
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Figure 5: A similar plot as figures |3] and H] highlighting the behaviour of all three 
extremal surfaces anchored at the same boundary points 9 = ±0.497r. Although two 
of the solutions (red lines) penetrate the dashed circle, these are not the minimal 
area and therefore are not part of the w{Da) set. 

The equation describing extremal surfaces can be found by minimizing the area 
of a spacelike codimension 2 surface given by 



leading to a second order differential equation for the function 6{r) which can be 
solved analytically in Z) = 3 dimensions or numerically for D > A dimensions. In 
any case, knowledge of 6{r) allow us to search for different solutions 6i{r) for which 
9i{oc) = 6*0, Vi, while rmm{Oi{r)) ^ rmin(6'j'(r)) for at least one i' . In other words, 
it allow us to look for extremal surfaces anchored at the same boundary points that 
do not share a common deepest bulk point rmin(6'j(r)) (figures H] and [5]) . 

The existence of such multiplicity of extremal surfaces with common boundary 
anchor points will in general preclude most of them from having minimal area and. 




) 



2 



(4.1) 
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as argued in |8J, can be used to show that no family of minimal area surfaces can 
cover the bulk in its entirety. Nevertheless, as a consistency check, we computed 
numerically the area for surfaces anchored at various boundary points and confirmed 
that, indeed, when more than one solution exists (with the same ^(oo)), the one with 
the bigger rmin has the smallest area. 

Therefore, our strategy to determine whether our charged boson star solutions 
display such hollow phases is to look for an extremal surface 6h{r) that obeys dh{oo) = 
7t/2 while having rmin(^/i) 7^ (figure H]) since, from the spacetime symmetry, we know 
that there exist one extremal surface for which ^(00) = 11/2 with r^[^{6) = 0. 

As a warming up exercise we fix the value of the parameters and q and, as we 
vary the central density of the scalar field ipo, "we observe the system to transition 
between solid and hollow phases as seen in figure |6l Note how on figure |6a] there 
is only one solution extending into the bulk and reaching a specific rmin for each 6, 
whereas on figure l6b] there is more than one value of rmin for a given boundary anchor 
point (a fixed ^0) near the r = region. 

{nn 
0.6 r 

0.5 - 

0.4 ; 

0.3 - 

0.2 ; 

0.1 - 

0.0 L 

0.0 0.2 0.4 0.6 0.8 1.0 0.0 0.2 0.4 0.6 0.8 1.0 

(a) V(0) = 0.2 (b) V(0) = 1.2 

Figure 6: A comparison of two different central densities for the plot of r^iT^{6{r)) 
in four dimensions, both cases with = and g = 0.1. It is clear how for anchor 
points roughly between O.lTvr < 9 <ti /2 there exist three distinct extremal surfaces 
with different values of rmin, however, only one of them has minimal area. See figures 
[31 m and E] for specific examples. 

Our goal for the next section is to determine the precise value iph of V'(O) fo^^ 
which this transition occur as a function of the parameters m? and q and, in the four 
dimensional case, compare it to tpc-, the central density threshold between stable and 
unstable configurations. 
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5 Phase Diagrams 



In this section we explore in further detail the conditions for which we observe solid 
and hollow phases. Since we are dealing with a one parameter family of solutions and 
have two free theory parameters, we should be able to construct a three dimensional 
phase diagram and find a two dimensional surface separating solid and hollow phases. 
To numerically accomplish this task we start by fixing q while varying m^. For each 
value of and q we look for the lowest value of V^(0) for which we observe multiple 
values of rinin(6'o) for the same, fixed, Oq and find a line separating the two regions, 
we then repeat this process for multiple different values of q. 

As seen in figures [7] and [HI we observe a similar behaviour for both three and 
four dimensional cases. Is is clear that the threshold value of decreases as we 
increase either or g, and a maximum, finite value is attained as the scalar mass 
approaches the BF bound and the charge goes to zero. 

1.0 h 




0.4 - 
0.2 - 

0.0 tiA I I I \ I I I \ I I I \ I I I \ I I I \ I I I \ 

-1.0 -0.8 -0.6 -0.4 -0.2 0.0 0.2 

Figure 7: Critical scalar field central density separating solid and hollow configu- 
rations for D = ?) dimensions with g = 0.1 in blue, q = 0.2 in light blue (dashed) 
and q = 0.3 in green (dotted). In all central density value below the line 

correspond to solid solutions, while above it lies the hollow regime. 

As discussed earlier in this paper, four dimensional boson stars are known to be 
unstable when the central scalar density rises above a critical value ipc in numerous 
different setups. Since we just established the existence of iph, the threshold between 
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Figure 8: Critical scalar field central density separating solid and hollow configu- 
rations for D = 4 dimensions with g = 0.1 in blue, q = 0.2 in light blue (dashed) 
and q = 0.3 in green (dotted). In all cases a central density value below the hne 
correspond to solid solutions, while above it lies the hollow regime. 

solid and hollow configurations, we want to compare it to ipc so we can determine 
whether the hollow solutions we find are in fact physically permitted. 

Stability of AdS charged boson stars and hollowed phases. 

The stability of boson stars has been a subject of intensive study in the past decades, 
while focus has been given to boson stars in flat spaces, similar results exist in AdS 
space. Despite subtle changes between the flat, curved, self interacting or charged 
cases, it is well known that four dimensional boson stars reach a maximum mass 
value for a finite central density ipc and are unstable past this point. The nature 
of the instability and how it depends on the various variations of the model, while 
important on their own, are not within the scope of this study, for us it suffices to 
know that the value of ip{0) = ipc for which the stars mass as a function of scalar 
central density, M{iIj{0)), is maximum represents the threshold between stable and 
unstable regimes and that this seems to be universal across different types of boson 
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star£l. 

In order to numerically determine ipc we fix q and and search for the highest 
value of the star mass M as a function of ipiO). Similarly to the phase diagrams above 
(figures [7] and [8]) , we find, for each value of q, a hue dividing stable and unstable 
regimes in the vs. i/jq plane. We observe that, within the range of parameters we 
studied, we always have ipc < i^h (as seen in figures [91 and [TOj) . indicating that the 
hollow solutions we found, while being perfectly fine in a mathematical sense, do not 
correspond to a physically preferred phase and suffer from dynamical instabilities. 

1.4 - 



1.2 - 




0.6 - 
0.4 - 
0.2 - 

0.0 L , \ , , , , \ , , , , \ , , , , \ , , , , \ , , , L 

-2.0 -1.5 -1.0 -0.5 0.0 

Figure 9: A phase diagram displaying both transitions found in D = 4 dimensions 
(stable — )■ unstable and solid — )■ hollow) with q = 0.1. The red line (below) is the 
stability threshold, a value of ipo above it (the light blue region) renders a dynamically 
unstable configuration. The blue line (above), once again, represents the transition 
between solid and hollow configurations. It is clear from this figure how, in the range 
studied, only solid configurations are physically allowed. 



-•^^For the interested reader we direct you the reviews [H] and [T^ and the work [TS^ for a lengthy 
discussion on boson star instabihties and more. 
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Figure 10: A phase diagram comparing the central densities V'c and V'/i as a function 
of for different values of q. The dotted lines are for q = 0.3, the dashed lines for 
q = 0.2, while the solid lines for q = 0.1. The warm coloured lines (below) correspond 
to transition between stable and unstable configurations, while the cold coloured lines 
(above) transition between solid and hollow phases. As we lower the charge q both 
ipc and iph increase, however their difference remains roughly unchanged, highlighting 
how hollow solutions are unstable for all range of parameters investigated. 

6 Final Comments 

In this paper we numerically investigated the behaviour of extremal, codimension 
2, spacelike surfaces in charged, asymptotically AdS, boson star backgrounds. Our 
main goal was to establish the conditions for which families of minimal area spacelike 
surfaces anchored on the boundary fail to fully cover the bulk of the spacetime. As 
discussed in the Introduction, this study was motivated by recent ideas regarding a 
possible connection between the holographic description of entaglement entropy and 
the gravity dual of a reduced density matrix as discussed in [8]. 

We observed that the relation between the star's mass as well the the star's charge 
and the central density of the scalar field in four dimensions behave much like what 
is known for both neutral boson stars in AdS and charged boson stars in flat space. 
Notably, the existence of a maximum mass for a finite ipi^) = V'c strongly hints 
towards the presence of a stability threshold and can be used to infer the physical 
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feasibility of the hollow solutions we were so interested in. In three dimensions we 
found results akin to what is known in the literature for other types of boson stars, 
in particular, we observed a behaviour similar as the one found in [16] for 2 + 1 
dimensional neutral boson stars in asymptotically AdS spacetime. 

Our analysis of the behaviour of extremal surfaces with fixed boundary points 
led us to the conclusion that, both in three and four dimensions, for fixed and 
q there is a maximum value for the central density of the scalar field iph for which 
the minimal area surfaces reach every point in the bulk space (figures [7] and [H]) • 
Therefore one should expect that charged boson stars with a high enough could 
provide a clear obstacle in the way of w{Da) being a universal candidate for R{A). 
However we saw that, at least in four dimensions, there is good evidence indicating 
that solutions with ^/'(O) > iph are unstable (figure [9]). Extrapolating well known 
results in the literature for both boson stars with and without gauge fields in fiat 
space, and boson stars without gauge field in AdS space, we find that for given w? 
and q there is a threshold central density value ipc for which the solutions cease to be 
stable if V'(O) > "^c- Remarkably, in the four dimension case in question, we found 
that for every pair of m? and g, ipc < ^h? i-e.: solutions for which w{Db) fail to 
cover the entire bulk and, in particular, z{Da) 3 w{Da) while w{Da) ^ z{Da), 
are physically unstable. Unfortunately, to the extent of this author's knowledge, 
much less in known about the stability of three dimensional boson stars, therefore 
precluding us from saying anything about the stability of both regimes we found. 

We believe the results found in this work support some of the ideas discussed in 
[S] and further explored in E] • The unstable character of hollowed solutions 

strengthens the proposal of w{Da) as a good candidate for R{A) and complements 
other recent works on the subject. We also believe that a deeper understanding 
of extremal surfaces on charged boson stars backgrounds can serve as a fruitful 
test ground for numerous holographic ideas including, but not restricted to, the 
holographic entanglement entropy, the holographic dual of a density matrix, zero 
temperature quantum phase transitions [T7|, etc. 
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